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Abstract 



The vacuum energy is calculated for a free, confornially-coupled scalar 
K* ' field on the orbifold space-time IR x S^/F where F is a finite subgroup of 

0(3) acting with fixed points. The energy vanishes when F is composed 

f~^ I of pure rotations but not otherwise. It is shown on general grounds that 

^D ' the same conclusion holds for all even-dimensional factored spheres and the 

^ . vacuum energies are given as generalised Bernoulli functions. The relevant 

Qv^ I (^-functions are analysed in some detail and several identities are incidentally 

derived. The general discussion is presented in terms of finite reflection 
groups. 
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1 Introduction 

Roughly speaking an orbifold, or a V-manifold |jl|] , is a quotient of a manifold by 
a group of symmetries with fixed points . The orbifold is singular at the images 
of the fixed points where, for example, the Riemannian curvature would diverge. 
Orbifolds can provide useful approximations to more complicated manifolds, such 
as Calabi-Yau spaces, occurring in string compactifications |^, ^, Q. They have 
also appeared in membrane compactification ||^. Some work has also been done 
on geometric quantisation on orbifolds 0. 

It is interesting to extend the usual theory of quantum fields on curved spaces 
to orbifolds. Some cases have already been considered. The ideal cosmic string is 
an example. We have earlier discussed some effects on the tetrahedron 0. 

In the present paper we consider a scalar field in the static space-time IRx S^F, 
where T is one of the finite subgroups of 0(3) associated with the regular solids. 
The fundamental domain of F is a curvilinear triangle or quadrilateral on S^ 
depending on whether F does or does not include reflections. The discussion is 
then extended to higher spheres. 

As usual, we consider free fields and, in the present work, as an example of 
a field theoretic calculation, we look at ways of evaluating the conformal scalar 
contribution to the total Casimir energy on the divided two-sphere. This will be 
sufficient to exhibit the various techniques and relations. 

Recently, there have been calculations of the vacuum energies on the hyperbolic 
space-times, IR x H^F, Zerbini et al 0. 

We have already shown that the total vacuum energy vanishes on the full 
sphere IRx S^ (and in fact on all even- dimensional spheres for conformal coupling). 
It will turn out that dividing by a purely rotational F does not alter this conclusion. 
Including reflections does produce a nonzero Casimir energy. 

The analysis is presented in some calculational detail because we think the 
discussion of the relevant heat-kernels and (^-functions has some general method- 
ological value. The calculation can be thought of as a rather explicit example of 
an equivariant situation, the general theory of which has been understood for a 
long time. 



2 Vacuum energy 

We consider an ultrastatic spacetime, IR x M, where the spatial section, M, is 
(i-dimensional. A scalar fleld satisfles the traditional equation of motion 

{n + CR + m^)if{x) = (1) 

with the D'Alembertian O = d^ — V^ . V^ is the usual Laplace-Beltrami operator 
on M. 



The general formula for the total vacuum, or Casimir, energy, E, for scalar 



fields in a static space-time is given in the form we need in [M. Thus 



^ = ^ lim L-'^'-'ka{s - 1/2) = ^ C.(-|) , (2) 

when this does not diverge, which it does not here. This will be discussed later. 
L is the scaling length. 

The ^-function, Q, is a Dirichlet series for the operator (— V^ + ^R + m^) on 
M; Cd{,s) = Y^x^~'^- Foi' constant scalar curvature, we can write 

C.(^) = f^ /"^ dr e-(-^+«^)^ t^-'K,{t) (3) 

1 [Sj JO 

where K^ is the integrated coincidence limit of the Laplacian heat-kernel on M. 
Incidentally, in the present case, the vacuum energy is also the negative of the 
effective (one-loop) Lagrangian, L^^\ which is related to the effective action, W^^\ 
by lyW =/dtL«. 

3 The Selberg trace formula on S^/F 

On the unit two-sphere, V^=— L^ where L is the usual angular momentum oper- 
ator. Then the defining equation for the heat kernel on S^ is 

{dr + L^) Ks2 (r, r'; r) = 6{r) 6{r, r') (4) 

which has the solution 

Ks2{r,r'-T) = ^e-'('+i)^Y<l>/^(r)Y<l>/*^(r') 

1,171 

1 OO 

= ^ E (2^ + 1) e-'('+^)" PKcos $(r, r')) 

where cos$(r, r') = r ■ r', r and r' being unit three-vectors. 

Next, the Mehler-Dirichlet integral representation for the Legendre polynomial 
is used to get 

Ks^ ir, V; r) = -^ £ (2/ + 1) e-'(^+^)^ f da /"^^+ ^^" , (5) 

2v27r^ ,^Q J* Vcos<l> — cosa 



which is quite standard e.g. ||TT| . 

We now restate our problem. The spatial manifold we want to examine initially 
is the two-sphere with a reduced symmetry, specifically M = S^F where F is a 
discrete subgroup of 0(3). F is either a cyclic (Cg), dihedral (D2q), tetrahedral 



(T), octahedral (O) or icosahedral group (Y). If one uses an extended group, T' 
(i.e. one that includes a reflection), then S^/F' is a fundamental spherical triangle in 
which the field satisfies Dirichlet or Neumann boundary conditions. For simplicity, 
we restrict the analysis initially to pure rotational F C SO (3). 

The method of images is used, with an argument similar to that found in 
McKean and Lax & Phillips |jl2|, |13|, to find the heat kernel on the divided sphere. 



The idea is to write the pre-image sum as a sum over conjugacy classes. It is easy 
to show that 

^s2/r(i") = / Y.^s'4r,-fr;T)dr 

= Ks2{r,r;T)dr + —Y,\{l}\Ks2{r,-fr;T)dr, (6) 

Jsyr |F| ^ Js2 

where 7 is a typical element, one of |{7}|, in the conjugacy class {7}. The elements 
of a class correspond to rotations through one fixed angle about a set of conjugate 
axes. For a given set of such axes, one corresponding class can be considered 
to be the primitive class, all other classes associated with these axes then being 
generated by this one. Thus the sum over all classes can be rewritten as a sum 
over primitive classes and the powers of these. Let q be the generic order of the 
rotation associated with the generic primitive class {7} so that 7"^ = id. Then 
I {7} I is just the number, Ug, of conjugate g-fold axes and we can write 



Ks2/r{T) = Ks2{r,r;T)dr + -—Y,nqY. Ks2{r,Yr;T) dr. 

Mr |F| Y ^1-^32 



(7) 



We can think of this rather obvious equation as a simple application of Selberg's 
trace formula to a compact, here spherical, domain. 

4 Evaluation of the trace formula 

We will now evaluate each term in the above trace formula, (^. The term corre- 
sponding to the identity element of the subgroup F, from (|]), is 

1 °° r r- 

^i:(2/ + l)e-'<«)7^^^^drj('dax,(2a). (8) 

where xK^^) = sin(2/ + 1)6/ sin 6 is the character of the /-representation of S0(3). 
Both integrals can be done to give 



1 

iTT 



^ (2/ + 1) e"'('+l)^ 



Up to scalings, this last quantity is nothing other than the extensively discussed 
rotational partition function of a gas of diatomic molecules |T^ , |Tl], |15[ . 

However we do not want this form for the first trace term (§), instead we shall 
leave it as follows 



2 



f;(2/ + l)e-'('+i)^/^d«xK2«). (9) 

7r|r|^ Jo 

Next we examine the second term in (^. A little geometry shows that, for 
7 = n(cu) eS0(3), 

cos $(r, 7r) = cos cu + (1 — cos cj) (r ■ n)^ . 

The summand of the double sum is then 

-| oo 



2v27r^ ,=0 



-1{1+1)t 



X I dr da sin(/ + ^)aRe , (10) 

^s^ ^0 ^c+(l-c)(r-n)2-cosa 

where c = cospu. We may choose n as the polar axis. 

Pulling the space integral through, to change the order of integration, and 
setting /i = cos 6, (|T0|) becomes 



^ V (2/ + 1) e-'('+^)" /" da sin(/ + i)a / dfi Re 
/27r ^ Jo Jo 



V^T^ 1=0 -^0 ^ -^0 Jc- cos a + (1 - c)fi'^ 



Since the radicand has to be positive, the two cases of c — cos a being either 
positive or negative must be considered. Define the quantity puJ = pcu ii pcu < n 
and puJ = 271 — puj if poo > tt. Hence 

1 v2 ^ ^^^ ^ ^^ -i{i+i)t 



7^ Vl - c ;=o 



^(2/ + l)e- 



{rpLxJ rl 1 

/ da sin(/ + i)a / , du , 
Jo J^^) J~^ 



Q{a) 

+ / da sin(/ + ^)a / d// _ =^ > 

where the terms in the braces can be evaluated to give 

< / da sin(/ + i)acosh$-l ('(5"^(a)) + /da sin(/ + i)asinh$-l (q 



2 (a 



The functions Q{a) and Q{a) are defined as Q{a) = (cos a — c)/(l — c) and 
Q{a) = (c — cosa)/(l — c). 

Integrating by parts leaves us with the desired, intermediate result, 



4^ ,/,,1^_„/■f - cos(2/ + 1)q; cosa 

cos 2a — cospuj 



71 f^^ Jo 



(11) 



The principal part, denoted by the P, is necessary in order to remove the 
singularity at a = puj/2 introduced by the partial integrations. 

Finally, summing over the cyclic group, u = 2ti /q, (excluding the identity, 
p = q), we get 



2f ,-««,.p/i,„-(2'+l)" 



vr 



1=0 



sma 



cot a — q cot qa 



(12) 



using a standard identity for the trigonometric sum. 
Noting the Jacobi theta-function relation. 



1=0 



I oo 

(2/ + l)a = J— Y. (-l)"e-("+"")'/^ 



and its derivative with respect to a, we find the heat kernel, from (|^), (|^) and 
(|l),tobe 



I I V / n=— oo 



I g-(a+7rn)7T 

da (a + nn) 

sin a 



r /"f g-(a+7rn)2/r 



sma 



cot a — g cot ga 



It is possible to subsume the outer sum into the integration range, producing the 
compact expression for the classical paths form of the trace formula. 



27]- Q'r/'i POO e""^/"^ 

Ks-^/T{r) = —- — -^P da — 

|i I (vrr)'^/^ JO sma 
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cot a — q cot qa 



(13) 



It should noted that this expression is akin to that obtained in the hyperbolic 



case by Donnelly []TB| and later by Balazs & Voros 

It is better to avoid the principal part by introducing a contour, C^, located 



just above, or just below, the real a axis and writing for (|T3D 



K~ 



Tx e 



t/A 



s2/r 



r 



|r| (7rr)3/2 



da 

Ca sm a 



a 



7 



Tin 



cot a — q cot qa 



(14) 



By expanding the sines as infinite sums of fractions and running Ca strictly 
along the real axis we can confirm that this integral is equal to the principal part 
form appearing in (|13|), [|18]. 

In the form (0) it is possible to effect a partial integration to get 

i^s7r(r) = p^X/«^|a(l-E-.)-^E-.^eotga|. (15) 
In the cyclic case, the term proportional to a goes away, leaving 

er/4 ^ g-a^/r 



2(7rr)^/^ JCc sma 



and (ITSI) can be written 



i^Q2/r(r) = — 7 



^s2/r 



r 



(17) 



Y.qn,K,{r)-{Y.n,-l)K,{T) 

_ 7 7 

where Ki{t) = Ka,2{r). 

There are several possible uses of these formulae, for example we could make 
an asymptotic expansion to give a power series in the limit r — i> 0, 

1 oo 

^svr(r) = - E Cut\ (18) 

''" fc=0,l,... 

This will be taken up later. The coefficients Ck are related to values of the C,- 
function given in the next section. 

5 The (^-function 

Important to the present work is a generalization of Riemann's zeta function. The 
incomplete Riemann (^-function (also known as Hurwitz's (^-function) is defined 
by 

oo -I 

Cb{z, a) = y^ — , for Re z > 1 . 

The analytic continuation of this, for all z, is 

, , , Til-z) r i-xY-^e-"'' , 

ZTTZ JCo I — e ^ 

where Co is the standard Riemann- Hankel contour around the positive real axis 



T9[| . In particular we will have cause to use the incomplete (^-function for the case 
a = |; so we give the contour representation 

c,(-).-!i^ni^/ (z^d.. ,19) 

1t:i J Co smhx 



Turning now to the field equation (||), the case of conformal couphng {m = 
0,^ = 1/8) is particularly simple. Since R = 2 for our scale, conformal coupling 
is equivalent to using the operator (L^ + 1/4) when finding the heat kernel and 
amounts to dropping the exp(r/4) factor in (|15D, which is the simplifying fact. 
Then the ^-function (j^) is the Mellin transform of the integrated heat kernel (jT^), 
up to this factor. We have 



Csvr(s) = 1^1 ^^. , / dtt^~' / da 



. g-a^i 



T\^/7TT{s) Jo Jca sin a 

X <a{l-J2'n'q)- ^J29^i^^^9^\ ■ (20) 

l 7 7 J 

Note that we have set t = 1/r. We will denote Cs2/r by Cr- 

The immediate objective is to invert the order of integration in ( pO]) so that 
the t-integral can be done. Some manoeuvring is necessary before this can be 
achieved because the t-integral will diverge when a is inside the 'light-cone' i.e. 
when x^ < y'^, where a = x+iy. We always require Re(— a^t) < for convergence. 

In terms of the real and imaginary parts, a = x + iy and t = u + iv, 

—a^t = — (x^ — y'^)u + 2xyv — i[{x'^ — y'^)v + 2xyu] . 

It is actually convenient to take the average of two contour integrals, one, Ca, 
above the real a-axis and the other, C'^, below. The q = I integral (corresponding 
to the full sphere) will then be 

1 „ f , ae " ^ 



Ciis) = dtt-^-' da- . (21) 

2^y^T\.{sjJo JCc+c'^ sma 

We study the behaviour as the a- and t-contours are rotated. 

Consider Ca rotated anti-clockwise onto the line x = y, then the t-integral 
would be along the negative imaginary axis of the complex t-plane. Moving Ca 
further round to the vertical, x = 0, would require the t-contour, Ct, to run from 
the origin to — oo below the negative real axis, which is a branch cut. Similarly, 
rotating C'a clockwise to the negative imaginary axis of a implies that the corre- 
sponding contour C'f- runs above the branch cut from t = to t = — oo. As the 
a-contours sweep passed the negative real axis, contributions are picked up from 
the poles but these cancel and so can be neglected. The next step is to combine the 
two t-contours, Ct and —Cj., to form one continuous path, C, that loops clockwise 
round the negative t-axis, passing through the origin. This is possible by reversing 
the rotated C'^ contour, which runs parallel, but opposite, to the rotated Ca one. 

These two rotated a-contours are next both moved onto the y-axis. Redefining 
a ^y ia and using the symmetry of the integrand in a we find. 



oo ^. ^a t 



r 1 r°° cy e 
Us) = -^^^J dtt-^-^ da^. (22) 

^^y^TL [s) Jc Jo smna 

The inversion of the order of the a- and t-integrations in ( p2|) is effected by 
deforming the t-contour, C, to run around the negative t-axis, avoiding the origin. 
Then we get 



Us) = -^w^ / d«-l— / dtt5-^e' 
iV^ri {s) Jo smha Jc 



This should be compared with (|2l|) . 
Integrating over t yields 



2r(|-3) roo C.2.-2 

Ci(s) = ^ , cosTTs / da 



^/7i^{s) Jo sinha 



where we must assume that Re s > 3/2 in order that the a-integration converge 
at the origin. 



Then, using (|l^), we finish with 



>.^ r(|-s)r(2s-l)cos7rs, , ,, ,, ,, 

Ci(^) = - ^' tIs)2^s-sJ^ Ck(2. -1,1) = 2U2s - 1, 1) . 

This is the conformal (^-function on the undivided sphere and is well known, follow- 
ing immediately from the eigenvalues, /(/ + 1) + 1/4 = {1 + 1/2)^ and degeneracies 
2(^ + 1/2), 

Cs2(s) = 2Cij(2s-l,i). (23) 

(See Dowker ["Vacuum energies on spheres and in cubes", unpublished 1983] and 
Camporesi [^ .) 

For the other terms in (pO[), it can similarly be shown that 

CO -^ POO C~'^ * 2^''~^ f^ Ci^^~^ 

dtt^'^^2 / da— cotga; = — — — — - / da ^ — cothga 



^/^TT{s) Jo Jo sin a r(2s) Jo sinha 

for Res > 3/2. Then the entire (^-function, (pO]) , can be written 

1 2^^^^ r^ o?^^"^ f 
^s^,M = Y\Yi2s)\o d«-i^|(2^-l)(l-E-.) + E^-.eothga 

' (24) 
valid for Res > 3/2. The same expression can be reached using the principal part 



form (|T3D 



In the case of the cyclic group, F = !^q, the first term in the braces in the 
integrand of (|2^) vanishes and the result simplifies to 



>2s-l ^nc-, Q/2*~^ 



Cs./.,(^) = 1X2^)70 ^"^i^^°'^^"- ^2') 



A simpler, eigenvalue derivation of this formula is given in the next section. 

Calling the ^-function of ( p5D Cq{^) for short, the general expression (^) reads 



(As) = ^ 



Y. l^gCgis) -(J2^i~ l)Ci(s) 



(26) 



where Ci(s) = (32 (s) is given by (|3|). 

For all r except the cyclic groups, J2^^q " 1 = |r|/2 and so we get 

Cr(s) = ^(s^s) + T^^Y.9ngCq{s) 

for r 7^ ^q. 

We now look for an analytic continuation of (P^. Expanding the hyperbolic 
cotangent in (pSj) in terms of exponential functions and integrating, yields pT| , 
p361 (3.552.1)] 

p=om=o (gm+p+ 2) 
Rearranging the sums, 

C,(.) = CRi2s, I) + -CRi2s -1,1)-^ E(2P + 1)Cr{2s, ^) , (27) 

1 1 p=0 

which would serve as an analytic continuation using the Hurwitz ^-function rela- 
tion ig |21|, pl073 (9.521.2)]. 

To produce a functional relation directly for the continuation of Cq{s), we 
consider the integral 



22^-ir(l-2s) 



/ da coth qa 

J Co sinh a 



2iTi J Co 

For Res > 3/2, this equals the integral in (|25|) and so is an appropriate continu- 
ation. 

For Res < 0, a large loop can be added without penalty to the contour Cq so 
forming a continuous path that ultimately surrounds all the poles on the imaginary 
axis. Then we get, from residues, and after rearranging the summations a little, 

^,... co...r(i-2.) ^ E -^(v)"' + P-i)i:(-ir(-")"" 

(mod q) 



Therefore, comparing with (0), the functional relation that must exist is 



2'-'r{2s){CR{2s,l 



1 "-^ 



y p=o 



E(2p+l)U2.,^ 



TC 



9-1 



q sm TTS 



«=i 



sm2£ 



E— ^ Ul-2.2)_22^^^(1_2, 



2g ' 



f28) 



which can be proved directly using the rational Hurwitz (^-functional relation and 
a couple of trigonometric sums found in Bromwich [^, p272 Exl8]. 
From (p4|), ( p7|) and (|28|) , there results the explicit forms 



(As) 



^Ci?(2s-l,|) + -^Egng 



y p=o 



2p+l- 



or, equivalently. 



(29) 
(30) 



22T(l-2s)cos7rs^ ^^ 1 



TT 



l-2s| 



P^'l--? 



a(l-2.,f)-22^C«(l-2.,^) 



2g 



Equations (|29|) and (|30|) also exhibit the analytic properties of the (^-function, 
Cp(s), which has only the one pole at s = 1 of residue l/|r| in accordance with 
general theory, which gives a residue of |S^/r|/47r. Further, the coefficients in the 
asymptotic expansion ( [T^ ) are given by the values of the (^-function at negative 
integers. The formula is 

U-k) = {-l)''k\Ck+i (31) 

It will be enough to give these values in the cyclic case. 



C.(-fc) 



1 



,2fc-i q-1 



qik + r 



B2k+2{l/2) + ^^—J2{2p+l)B2k+i{{2p+l)/2q) (32) 



p=0 



in terms of Bernoulli polynomials. 

The important value Cr(0) is given in general by 

Cr(0) =Cd/2 -no 



(33) 



where no is the number of zero modes. For conformal coupling uq vanishes. From 
(0),or (^, we find {d = 2) 



CM 



2q^-l 
12q 
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We also note, most easily from (pOD, that (^.{s) has zeros at s = —(2k + l)/2, 
with fc = 0, 1, . . ., in particular at s = —1/2. Thus, from @, the vacuum energy 
vanishes, as mentioned earlier. It also vanishes on space-times with the structure 
IR X IR^'^x S^/r. Of course, it must not be forgotten that this conclusion depends 
upon the assumption of conformal coupling. 

Expressions for the coefficients Ck can be found by expansion of the trigono- 
metric functions in (|TB|) or, more easily, in (|TB|) in powers of a. The result appears 
as a polynomial in q whose coefficients are products of Bernoulli numbers. Equiva- 
lence with (|32D devolves upon an interesting identity and is discussed in Appendix 
A. 



6 A direct eigenvalue method 

In this section we describe the more direct eigenvalue method of deriving the 
C-functionand will later enlarge the discussion to include the extended groups. 
These are easily allowed for by reinterpreting F in (^ to be an extended group and 
inserting phase (sign) factors in the sum such that the direct elements, i.e. those 
with an even number of reflections, have a factor of +1 while those elements with 
an odd number of reflections have a factor of -|-1 or —1 depending, respectively, 
on whether we require Neumann or Dirichlet boundary conditions at the edges of 
the fundamental domain. 

On Syr, the eigenvalues of (L^ + l/A) are ti^/A for n = 21 + 1 and / = 0, 1, 2, . . .. 
Iid{l), = d{n), is the degeneracy of each eigenvalue labelled by I then, by deflnition, 
the conformal (^-functionon S^/F is 

C(,s) = 2- E ^. (34) 

n=l,3,... "■ 

From flrst principles d{l) is the number of times the trivial representation of F 
occurs in the /-representation of SO (3), (c/ p3| ) 



1 2/ + 1 1 '^""^ 

^(0 = TfrX^,il) = -T^ + TfrX''^T.xXr) (35) 

|J- I 7 |J- I |J- I 7 p=l 



where the cyclic character is 



sin((a + l).p/,) 

sm(7rj9/g) 



Thus for F = !^g (which is sufficient) 



l.-i 



^.(0 = -Ex,(7^), (37) 

y p=o 
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which reduces to 2[//g] + 1, where [//g] is the integer part of l/q. 
Evaluation of the (^-functionfrom (|3^ proceeds as follows. 



2^' 1^ I sin [mip/q) 



^'^'^ l^oM-pMn^s,. 



n 



2s 



= ^ V 1 V ( ^^^(^^P/l) _ sin {2mxp/q) \ 

q ^, sm{7rp/q) ^J^,^ \ n^^ (2nY^ j' ^ > 

Separating off the p = term and then setting p —>■ q — p we find 

> . s 2 , ,^22"^:^ 1 ^ sin (mrp/q) , , 

q g ^1 sin(7rp/g) ^^2,... ^ 

Next we recall Hurwitz's formula 

.^-^ sin(2ra7ra) (27r)^'' , , , , xx / x 

to rewrite (^) in terms of Ci?-f unctions, 
Cg{s) = -Ch{2s- 1,1/2) 

+ A vfo ^ ■ ^ • , / , C/?(l - 2s, f) - Cn{l - 2s, ^) . (41) 

4gr(2s) smvrs rr^ sm(7r]9/g) L ^^ "^'^ -i 

We could stop here but, if agreement with our previous result, (|30|), is required, 
the relation, 

can be used. 

It is easy to introduce a U(l) twisting into the character formula, (0), by 
writing 

d,{l, t) = - E ^^^((2^ + ^)^P/g) cos(27rpt/g) = rf,(/ + t) + d,{l - t) , (42) 
Q. p=o sm(^7rp/gj 

where t is an integer between and q. One can think of this as the effect of an 
Aharonov-Bohm flux along the rotation axis. The modification will clearly carry 



through unchanged into the final formula, like (H), for the (^-functionwhich as a 
consequence will have the same zeros as before. The same conclusion holds for 
the twisting, 

^(0 = p^Ex.(7)x;(7) = ^^^^^^T^ + p^E^/Ex.(r)x;(7^), (43) 

|J- I 7 1^ I |J- I 7 p=i 
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that results from taking the field to belong to any representation, D, of F. 

It must be said that there is no substantial difference between the method of 
this section and that given in section 5. Not surprisingly, the integral ([TT| ) can 
be done to yield the expression ( P?] ) for dq{l). We have taken the route we have 



in order to exhibit various aspects of the analysis and to provide checks of the 
formulae. 

It is interesting therefore to make contact with the work of Polya and Meyer 



2% who gave Molien generating functions for the degeneracies of each harmonic 
basis of degree / invariant under a point group of three-dimensional Euclidean 
space. (See also Laporte pSJ .) 



The generating function is defined to be 



h{a)=Y.d{lW. (44) 



From (|37|) , a geometric summation produces the cyclic form 



haia) 



:i+a)^' 1 






'' q f^ 1 - 2(T cos{2Tcp/q) + cr^ 



showing the connection with the Molien determinant expression p4|, A further 
summation yields, c/ [^ Table 1, 

1 — cr 1 — cr'J 

Expanding gives again dq{l) = 2[l/q] + 1. 

It is amusing to note that the generating function is, to a factor, the integrated 
heat kernel of the (pseudo-) operator (L^ + 1/4)^/^, 



K^/'^{t) = Tr exp ( - rJL^ + 1/4) 



This can be seen by setting cr = e "^ so that, making redefinitions where 
appropriate, 

h{T) = ^rf(0e-"' = e"/2 ^ rf(n)e-™/2 

/ 71=1,3,... 



e'/^Tr exp(-rVL2 + l/4) 



where the trace is over all eigenstates on S^/F. 
The C-functionsare related by 

C(.) = C^(.) = Tr A-^ = Tr (A'/^'-^ = C^^ (2s) = C'^\2s), (46) 
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i-is) Jo 



with 

ry^(.,) = 

T{s) Jo 

where, by the above, 

K'/\T) = e-^/^h{T). 

For the cychc case using ( 145|) 

/ cothgr^^ 

" ^ ^ 2sinhr/2 



(47) 



Therefore, for Re s > 3, setting t = 2a 



2*~^ r°° fy'^~'^ 
Cl/'{s) = —- da— — cothga, (48) 

1 (s) Jo smh a 

and, using (^6]), we regain the answer given above, ([25|), as we must. 

The general result (^^ or (^6D also follows using (^) and also is in accordance 
with the equation on p. 139 of [^. The ^-functions for the finite subgroups of 
0(3) thus take the compact integral form, written using the Hankel contour so as 
to hold for all s, 

(As) = ^^^^^^X/-(-)""^y'M' (49) 

where the K^^'^{t) are given by 

^1/2. ^ ^ cosh(dor/2) 

^ ^ ' 2sinh(rfir/2)sinh(d2r/2)' ^ ^ 

with the constants, (^0,(^1,^2) equal to (g, g, 1) for Wig, (2g + l,2g, 2) for D2q, 
(6, 4, 3) for T, (9, 6, 4) for O and (15, 10, 6) for Y. 

We have earlier shown that the Cj^{s) vanish at s = —{2k + l)/2. A; = 0, 1, . . .. 
This can also be seen very simply as follows. 

At s = —{2k + l)/2, the integrand of (|49|) is single-valued and the contour can 
be shrunk down to a small loop around the origin reducing the computation to 
one of residues. The evenness of the integrand, (0), means that the residue at 
the origin is zero without the need for any particular representation. 

In contrast, the C-functionsfor the eaiended groups do not vanish at s = —{2k+ 
l)/2. In the case IR x S^/F' we have explicitly 

^ 4smh(air/2) smh(a2T/2) 

with the corresponding constants, (rfi, ^2, c^o); and find for the vacuum energy. 
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The sign depends on whether one chooses Dirichlet or Neumann conditions, which 
itself is determined by setting the twisting, x(^)) equal either to det A or to 1, as 
described earlier. 

The specific numerical values are ±11/192 for T', ±29/256 for O' and ±89/384 
for Y'. (See Table 1.) 

7 Higher spheres and reflection groups 

We enlarge the analysis to space-times of the form IR x IR^ x S^'/F and also give 
a more systematic account of the extended groups. 
In place of (H) there is now the general equation 

E = , )^, ; rrXdk-k - 1/2), (53) 

2(27r)'=(2A;±l)!!^^ ' '' ^ ' 

when finite, for the vacuum energy per unit volume in the IR^ space. 

We assume that, for all k, the scalar field is conformally coupled in the static 
space-time, IR x S'^/F, so that the relevant (conformal) eigenvalues are still (n ± 
d-2)/2,n = 1,3,.... 

In the present case (d is Cr given by (H^) where the heat kernel is related to 
the generating function by 

irV2(^)^g-{d-iW2^(^)_ (54) 

From (111) it follows that 

2fc+l 

r=0 



C(-A:-l/2) = -(^] e-('^-^W^Mr) 



Hence the behaviour of the combination e~*^'^~^^'^/^/i(r) = a^'^^^^^'^h(a) under r -^ 
—T, i.e. under a -^ ^/cr, is relevant; if even, the residue is zero. 

The structure of the (twisted) harmonic generating function, h^{a), on factored 
spheres, S7r, is given by, 0, |6|, 0, 

ft^M=(l^^_4lfi_, (55) 

^^ ' |F| "^ det{l-aA)' ^ ' 

which incorporates Molien's theorem. The situation now is that F is a finite 
subgroup (extended or not, as the case may be) of 0{d + 1) and the A are its 
{d + 1) X {d+1) matrix representatives. A twisting, x{^)i has been included, x(^) 
is shorthand for x{l) where A represents 7. Further, since A is real orthogonal, 
x(v4) could be replaced by its real part. 
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hx{cy) is a rational function but, at the moment, we do not need to be more 



precise than (|55|) . Then it is easily shown that 



^(d-i)/2(i _ ^2) ^ ^.(^)(i _ (_i)^det A) 

|r| ^ det(l-aA) • ^^^' 



For even-dimensional spheres, the sum in ( pq) reduces to one over those ele- 
ments of r that contain an odd number of reflections so immediately we see that 
cr(a'-i)/2/^^^fj^ is an even function when F contains only rotations. This shows quite 
generally that the (^-functionson such factored, even- dimensional spheres vanish 
at negative odd integers, (for conformal coupling of course). For odd spheres only 
the pure rotations contribute to (|56|). 

The generating function (^) is an example of a Poincare (or Hilbert) series 
and its derivation is an episode in Invariant Theory ||28|. For a finite reflection 



group, which is another name for an extended group, the structure of the usual 
harmonic generating function is just (pSj) with x(^) equal to unity. 



M^)-^i:5^j(ibi)-(^--^)n(r^. (57) 

where the di, i = 1,2, ... ,d + 1, are the degrees of the algebraically independent 
generating members of the invariant polynomial ring associated with the action 
of the group, F, on the d+1 dimensional vector space. Without loss of generality, 
F may be taken to be a subgroup oi 0{d+ 1). 

Since there is always the basic invariant r'^,= x"^ + y"^ + z"^ + ■ ■ ■, of degree 
2,= dd+i, we see that the corresponding term, i = d + 1, in the denominator 
of ( p7[ ) is cancelled by the 'harmonic factor', (1 — a^), leaving just the degrees 
(di, d2, . . . , dd) = d. 

The properties of reflection groups were worked out by Coxeter |^ , Chevalley 
30[| , Shephard and Todd [^], Steinberg [^ and Solomon |0. Flatto gives a read- 



able review in |^. An introductory account can be found in Grove and Benson 



35[]. For a more advanced expository treatment see Humphreys [B6|. Physicists 



might care to consult [^. Carter, |3^, contains useful algebraic information. 



Some standard, relevant facts follow. 

For small d, the degrees di can be determined by using an identity relating 
them to the geometrical properties of the group action. For example, it can be 
shown that Hi di is the order of the group and that J2i{di — 1) equals the number 
of reflections. In the two-sphere case these facts rapidly lead to the numbers cited 
earlier and we see that do = di + d2 — I is the number of reflections. 

Rewriting ( p7| ) in terms of the time r we find 

fJl 2 smh(dir/2) 
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Note that the effect of the conformal couphng term, e'^''"^''^/^, is here associated 
with the number of reffection hypersurfaces, the number of 'positive roots'. Using 
(0) we can check that (|58D reproduces the heat kernel (|5TD on S^/F'. 

On restricting the reflection group to its rotational subgroup, another gen- 
erating invariant polynomial comes into play. This is the Jacobian, of degree 
do = J2i=i{di — 1), of the original generating set. Under the action of F, J 
changes to {detA)J. It is an odd form and, as such, is the product of do linear 
factors, the vanishing of each of which is the equation of a reflecting hyperplane. 

Multiplying invariant polynomials by J produces 5A;ew>invariant polynomials 
and conversely all such polynomials have this product form. Thus their generating 
function is 

^*-nT^-^i:5^J^. (59) 






1=1"- " r I A 

This can be referred to as the Dirichlet generating function since the odd elements 
of F enter with a minus sign, cf |^^ . The harmonic generating function is 

h,{r) = e(^-W2e-^o-/^ {[ . } (60) 

fj^ 2 smh{diT/2) 

The corresponding Neumann function has already been given in (^). 

Since skew-invariants are strictly invariant under pure rotations, the generating 
function for F C SO{d + 1) is 

(1 + "'•) n T^ = e''^'>"' ~*(*-/2) n j-^y («) 

Multiplying by the harmonic factor, (1 — a^) = 2e~^ sinh r, removes the {d + l)-th 
term from the product and yields 

Mr) = kjr) + ft„(r) = 2e<-')'/= cosh(4r/2) Q j-j^. (62) 



From (^) and (|H^) we can check (|SD| ) for the heat kernel on S^/F in the 
rotational, SO (3) case. 

Equation (^) with ( pT| ) is a particular case of Barnes' C-function[|^ encoun- 
tered in our earlier work |^. We can write in the S^ case 

Cp(s) = C2 (2s, (rfi + d2- do)/2 I rfi, da) = C2 {2s, 1/2 | di, ^2) . (63) 

The general definition of Cd(s, * | *) is 

^ ( \A\ ^r(l - s) f exp(-ar)(-r)^-^ 
Cd{s,a\d) = ^^ Ldr—^ 



2vr Jl nti(l-exp(-ci.r)) 

00 -1 

Y- — , Res>d (64) 

^o(a + m.d)^' ^ ' 
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where d = (di,d2, ■ ■ ■ , dd) and m = {mi, m2, . . . , rrid) with the rrii non-negative 
integers. The location of the infinite contour L depends on the relative positions 



of the numbers a and d, [^ . In the present case it can be taken to be the Hankel 
contour, Cq. 

Barnes gives the values of Cd(s, * | *) at all integral values of s. In particular 
the values at the negative integers are the generalised Bernoulli functions which 
result from a residue evaluation via an expansion of the integrand. In particular 



in the notation of [^. Here, and in the following, the sums and products run 
from i = 1 to i = d. 

The Bernoulli polynomials are developed most thoroughly by Norlund El]. As 



noted by Hirzebruch |^ they are essentially the same as the Todd polynomials. 



In the general case, constructing the heat kernel one gets 



in) ^^2sinh(dir/2)' ^ ^ 



According to (|5^) with (57), we find that the Neumann and Dirichlet (- 
functionson S^'/F are given by (Q). Specifically, 



C(s) = Cd(2s,(d-l)/2|d), (67) 

C(s) = Cd {2s,Ed^ -{d- l)/2 I d) , (68) 

generalising the previous expression (^). 

We already know that the Neumann and Dirichlet vacuum energies on IR x 
IR^ X S'^/r are opposite in sign for even spheres and of the same sign for odd 
ones. Equivalently, we could use 

5f)(Erf. - a; I d) = (-l)"5(™)(x I d) (69) 

and the answer is 

(—1 \ku] 

where we have used 2Y{di = |r|. Here, and later, the upper sign is for Neumann 
and the lower for Dirichlet conditions and we recall that the field is conformally 
coupled in only (1 + (i)-dimensions whatever the value of k. 
The vacuum energy on IR x IR^ x S^/F is 



, , (-l)^fe! p(2) 



-^q2/r = ='= _t./oL I oM -^2fc+3 (1/2 I di, d2) ■ 



For k = 0, using 



i?f(7|a,/5) 
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(27 - a - p) (27(7 - a - /5) + a/3) 



(p^) can be checked. 

The rather trivial case of the one-sphere, i.e. the space-time H x S^/[g], gives 
the vacuum energy E = — g/24 for both Dirichlet and Neumann conditions. For 
the circle (g = 1), adding the two values produces the standard periodic result, 
E= -1/12. 

For the regular tessellations of the three-sphere, we find the values in Table 1. 



Group 

r 


Degrees 
d 


Casimir energy^ 

E 


[5] 


(q) 


-g/24 


[3,3] 
[3,4] 
[3,5] 


(3,4) 
(4,6) 
(6,10) 


^11/192 
^29/256 

^89/384 


[3,3,3] 
[3,3,4] 
[3,4,3] 
[3,3,5] 


(3,4,5) 
(4,6,8) 
(6,8,12) 
(12,20,30) 


-601/28800 
-3557/46080 
-69391/414720 
-3178447/5184000 



"Where there are two signs, the upper is for 
Neumann and the lower for Dirichlet boundary 
conditions. 



Table 1: Vacuum energies for the reflection groups on S^, S"^ and S^ 



8 The heat kernel expansion 



In place of ([T8|) we have 



K 



S'^/r 



J) 



r 



d/2 



E c,t' 



(71) 



fc=0,l/2,... 

where any half odd-integral powers are due to the boundary conditions. 



19 



Equation (^) becomes 

C,(-A;) = {-lfk\Cu+,/2 (72) 

and enables us to find Ck+d/2, k an integer. A further result we need is the 
following. 

On a d- dimensional manifold with boundary, Cd{s) has poles at s = (d — m)/2, 
for m = 0, 1, . . . , d — 1 and m = d+ l,d + 3, . . ., with residues 

r((d-m)/2)^"^/^- ^^^^ 

The possibly nonzero heat kernel coefficients are thus Co,Ci/2,Ci, . . . ,Cd/2 
and Cd/2+i, Cd/2+2, ■ ■ ■■ 

The importance of the coefficients, so far as field theory goes, is that in a space- 
time of the form IR x IR^ x Md, C(i+2fc+d)/2 is proportional to the coefficient of the 
pole term in the effective Lagrangian and a nonzero value indicates an ultraviolet 
divergence. Our interest here lies in even 2k (most typically zero), since, even 
if no divergences appear, we do not wish to undertake the technical evaluation 
of the derivative of the C-functionthat necessarily appears in the vacuum energy. 
(The result can be expressed as a multiple F-function, [^ .) 

The issue of ultraviolet divergences thus amounts to the existence of a pole at 
s = —k — 1/2, as intimated in section 2. The Barnes C-function(0), Cd{,s,a \ d), 
has poles only at positive s,= 1,2, . . . ,d — 1, so there is no ultraviolet divergence 
and no conformal anomaly. Expressions (^ and (^) are thus sensibly ffnite, 
although, de facto, evaluation is really enough, as we have seen. 

We calculate the coefficients Cd/2+k, k = 1,2, . . ., using ( |7^ ) and ([651). From 
(^), the residue this time involves the even part of K^^'^(t) and we need 

a(^"i)/2/^^(a) + a(i-^)/2/i^(a-i) = 

^(^-i)/2(i_^2) ^*(^A)il + i-irdetA) 

|F| Y det(l-aA) ' ^ ' 

If d is even, only pure rotations elements survive while, for odd d, only odd parity 
ones remain; in particular the identity does not contribute for odd d. 

From (^) and (^) expressions for the coefficients in terms of generalised 
Bernoulli numbers can be obtained 

^'/'^' = ^^^' |rp+Tfc)!fc! ^^+^^^^^ " ^^/^ ' '^^' (^ = 1' 2' • • •)• (75) 



The coefficient Cd/2 is given by (^ 



Cd/2 = {Tr^Bi^\{d-l)/2\d). (76) 
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The remaining coefficients can be found from the residues of the Barnes (- 
function 

Ci,-m = (T)^'^-'^ |r|(/_^[y(^j_^), ^g.((^ - l)/2 Id), (A: = 1, 2, . . . , rf). 

(77) 
For odd- dimensional spheres, and purely rotational F, we thus generalise the 
result that the heat kernel asymptotic expansion on the full sphere terminates, in 
this case with the Cd/2-1/2 term. 

A particularly interesting example is the three-sphere since we can make con- 
tact with the results of Ray |^ on the relevant generating functions. 

The asymptotic expansion of the heat kernel on orbifolds has been discussed 



by Donnelly [|T^, ^, and Briining & Heintze p6[ amongst others. The particular 
case of factored spheres is considered in some detail in |^, the structure of the 
Poincare series being presented more abstractly than in [2^ . 

The fact that the eigenvalues are proportional to squares of integers means 
that ( |i7| ) can be used immediately to give a series for the free energy, F, of the 
finite temperature field theory. Thus, in general, 

where (5 = 1/kT. 

In the cyclic case, on S^, 

F(3) = J-Y- coth(gm/3/2) 
'^^ ' 2(3^^msmh{m(3/2y 

9 Other couplings 

It is also possible to perform the calculation for massive or for minimally cou- 
pled fields but then Bessel functions would arise and the analysis would be more 
complicated(c/[|,0). 

The minimally coupled case is important as it corresponds to conformal cou- 
pling on S^/F regarded as a Euclideanised space-time. Interest would then centre 
on the evaluation of C' (0) s-nd of Cr(0)5 which is proportional to the conformal 
anomaly. In fact it is easy to evaluate this last quantity using the results of the 
present paper since there is no extra difficulty in finding the heat kernel expansion 
coefficients. The relation between the conformal heat kernel and that. Km, for 
the general equation of motion (|I]) is 

Kmir) = expr(m2 + ^((4^ - l)d + l))K{r) (79) 
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which, upon expansion, allows the new coefficients to be determined in terms of 
the old in a well known way. 

In the case of two dimensions, and minimal coupling, we find the conformal 
anomaly 

C(0) = (Ci - -C^) = — i— (12(rfi + d^f - Adid2 + 60(rfi + d2) - 45). (80) 
4 6bO\l I 

The space-time region can be thought of as a curvilinear triangle, and the 
conformal anomaly can also be found by conformal transformation methods. 



10 Conclusion 

We have found that, if F contains only rotations, the total vacuum energy, for 
conformal coupling on the divided even-dimensional sphere, S^'^/F, is the same as 
that on whole sphere, namely zero. One would expect, however, that the vacuum 
average of the /oca/ energy density, (Too (a:)), would be nonzero and would diverge 
as the fixed points are approached, just as for a cosmic string. 

If F contains reflections, i.e. is an extended group, then the vacuum energy is 
nonzero. In this case, because of the Dirichlet or Neumann boundary conditions, 
we would expect the local density to diverge at the sides of the fundamental 
triangle. 

The vacuum energy on an odd-dimensional sphere is always nonzero, since, 
this time, it is the pure rotation elements of F that contribute, and these must be 
present. The curiosity here is that Neumann and Dirichlet boundary conditions 
give the same result, just half that of the strictly rotational F C SO{d + 1) case, 
the half factor being simply a volume effect. 

A number of extensions suggest themselves. It would be straightforward to 
evaluate the finite temperature corrections and also to include the effect of a 
monopole magnetic field. In the latter case, by adding, in an ad hoc fashion, an 
extra term proportional to the square of the magnetic field to the Hamiltonian, it 
is possible to make the eigenvalues perfect squares so that the generating function 
can again be used. The only difference is that the angular momentum quantum 
number starts at the monopole number. 

A further extension that might be considered is that to the Dirac equation. 
However there appears to be an obstruction to setting up spin-half theory on S^/F, 
at least using images. 
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A A Bernoulli identity 

In the limit as r — > 0, the heat kernel can be expanded as a power series: 



1 oo 
'^ k=0 



We will just concentrate on the cyclic case, T = 'Eg. 

As mentioned in section 5, the trigonometric functions in the integral repre- 
sentation of the heat kernel, ([T3|) or ([iBD are written as power series in a involving 
Bernoulli numbers, and the integration then carried out. One finds, cf [|^, 



qi (^ „ II. 

^ U ^ ^ \2n)\{2m)\ (n + m + 1)! ^' ^ ^ " ^ 

Manipulation of the double sum yields the heat kernel coefficients, for A; > 0, 

^^ ^ qk\ |"^2fe(|) + ^^— Y Y. (Im) (^^™ - ^)B2k-2m{\) B^mj ■ 



These coefficients give the values of the zeta function at negative integers, via (|3T|). 



The Ck can also be found directly using the values of the Hurwitz zeta function, 
CR{—k, b) = —Bk+i{b)/{k + 1) for k >0, and the result is given in (|32D. 

The relation needed to reconcile these two expressions is the Bernoulli sum 
identity, 

- 2p+l ' 



E Q (9"" - ^)B2k-2n.C^) B^^ = k q''-'Y.i^P + 1) 52fc- 
m=0 p=0 



H 2g 



It is simple to verify this relation for k = 0,1 while for A; > 1 we have recourse to 
Lemma 1 of Apostol, [^, which says that for r > 1, there is the identity. 



1/-I 2r 

2rJ2^ B2r-i{t) = E f;) B, B2r-t y'-' + (2r - 1) u^-^^B^,. (82) 



This can be used to get 

£ (2p + 1) B„r-^) = 4 E (,^^) 52™ 52.-2,n q'~^-{2'-^- - 2) 

p=0 ^'^ m=0 

2k -1 



- ^2-2k(o2-2k 9\ p 

H ^7— 9 \2 - Z) B2k 
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Thus 



E 

p=0 



kq''^-'Y.(2p+l)B,,_,C-^) = {2k-l)B,,{l)+Y: © q"'-'"'B,,.,^B,Ul) , 



k 
m=0 



and we are left with having to show that 



k 2k 

{2k - l)B2k{l) = - E (2m) ^2fc-2m(|) B2m = E (m) (^ ~ 2^"™) B2k-mB„ 



m=0 



771=0 



Again the identity ( ]82D from Apostol, [^, gives 

/2A; 2A; — 1 ^ 

i2kB2k-i - {2k - l)B2k) - 2 i—B2k-iil) - -l^2^B2k 

for the right hand side of (|83D . This simphfies to 

{2k-l){2'-''-l)B2u 
and so we have proved the BernouUi sum, (pi]). 



(83) 
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